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ABSTRACT 


This  investigation  uses  ballistic  theory  in  the  analysis  of  the  behav¬ 
ior  of  an  electron  beam  passing  through  alternate  gap  and  drift  regions, 
with  the  gap  regions  having  excitation  fields.  The  analysis  develops  the 
response  to  a  complicated  frequency  spectrum  of  the  drive  signal.  First 
the  double  frequency  case  is  analyzed.  Relativistic  effects  are  taken  into 
account  and  their  influence  on  the  current  response  studied.  Then  the 
theory  is  extended  to  the  more  complicated  case  of  a  Gaussian  spectrum. 
The  first-order  bunching  theory  is  used  to  plot  current  response  curves. 

An  estimate  of  the  pulse  distortion  resulting  from  nonlinear  electron  beam 
dynamics  is  obtained  from  the  curves.  It  is  also  of  interest  that  the  envel¬ 
ope  shape  of  the  exit  current  is  almost  completely  independent  of  the  r-f 
frequency.  The  large-signal,  finite-gap  analysis  is  carried  out,  and  the 
results  extended  to  the  multiple -cavity  klystron. 


I.  INTRODUCTION 


A.  OBJECTIVE 

Many  stages  of  development  followed  the  invention  of  the  klystron 
by  the  Varian  brothers  in  1939.  Until  now,  most  of  the  analysis  and  dis¬ 
cussion  has  been  limited  to  the  case  of  a  continuous -wave  drive  signal. 

The  purpose  here  is  to  open  the  door  to  investigation  of  the  numerous  prob¬ 
lems  associated  with  the  use  of  pulsed  microwave  amplifiers  for  amplifi¬ 
cation  of  nanosecond  pulses  whose  pulse  lengths  are  of  the  order  of  several 
cycles  of  the  r-f  carrier  frequency.  The  main  concern  is  with  high -power 
amplifiers  with  average  power  capabilities  comparable  to  conventional 
pulsed  amplifiers;  thus  for  comparable  repetition  rates,  the  peak  power 
would  be  higher  by  the  inverse  ratio  of  the  pulse  lengths. 

Several  factors  may  be  important  in  determining  the  pulse  response 
capability  of  high-power  amplifiers- such  as  the  klystron.  It  follows,  from 
Fourier  analysis,  that  a  long  pulse  of  constant  carrier  frequency  includes 
a  narrow  bandwidth,  while  a  pulse  that  is  short  in  terms  of  cycles  of  the 
r-f  carrier  has  a  broad  frequency  spectrum.  The  spectrum  of  the  long 

signal  can,  however,  be  significantly  broadened  by  introducing  modulation. 

3 

Klauder  showed  that  to  utilize  the  transmitting  tubes  efficiently,  this  modu¬ 
lation  must  take  the  form  of  frequency  modulation.  By  this  method  one  can 
introduce  the  frequency-spread  characteristic  of  a  short  pulse  within  the 

4 

envelope  of  a  long-duration  signal.  Klauder  also  showed  certain  advantages 
of  short  constant-frequency  r-f  signals  over  the  long  signals  with  linear 
frequency  modulation.  This  emphasizes  the  importance  of  nanosecond  pulse 
studies.  One  of  the  problems  that  arises  is  that  the  broad  frequency  spec- 


-1- 


trum  associated  with  the  short  pulses  might  be  affected  by  the  bandwidth 
of  the  circuits  associated  with  the  amplifier,  which  will  limit  the  response 
and  therefore  cause  distortion  of  the  pulse.  Operation  of  the  amplifier  at 
maximum  efficiency  entails  nonlinear  behavior  in  the  electron  beam  dy¬ 
namics.  This  will  cause  the  frequency  spectrum  of  the  output  pulse  to  be 
altered  from  that  of  the  input  pulse  and  produce  distortion.  Ballistic  theory 
will  be  utilized  to  determine  the  response  of  klystrons  to  the  complicated 
frequency  spectrum  and  lead  to  an  estimate  of  the  pulse  distortion. 

Studies  in  this  direction  will  provide  a  solution  for  the  conflicting 
requirements  of  long  range  and  high  resolution  in  radar  systems.  Re¬ 
solution  depends  on  the  transmitted  pulse  bandwidth,  and  nanosecond  pulses 
will,  no  doubt,  satisfy  the  conditions  for  high  resolution.  For  long-range 
capabilities,  large  power  requirements  are  necessary;  hence,  high-power 
nanosecond  pulses  are  expected  to  solve  the  two  conflicting  radar  require¬ 
ments.  Radar  systems  that  yield  simultaneous  information  about  the  range 

4 

and  velocity  of  a  target  would  be  useful  in  certain  applications.  Klauder 
showed  an  inherent  ambiguity  in  a  simultaneous  determination  of  both  the 
range  and  velocity  of  a  moving  target,  when  using  the  so-called  "chirp" 
scheme.  If  the  transmitted  signal  with  an  ambiguity  function  that  is  highly 
peaked  only  at  about  t  =  0  exists,  then  high  resolution  is  expected  in  both 
range  and  velocity.  Using  analogs  from  quantum  mechanics,  Klauder  showed 
that  the  sequence  of  signals,  f(t),  that  satisfy  these  conditions  are: 

1/4  -  4- 

f(t)  =  ,  ^  ■■  -  H  (n/^)  e  2 

J  1/2  ,  ^n  ^ 

qtr  '  nJ  2 
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where  H^(z)  represents  the  Hermite  polynomial  defined  by 

H  (z)  =  (-1)’^  e^^  il.  e-=^^  . 

“  j  n 

dz“ 

When  n  is  taken  equal  to  zero,  the  Gaussian  envelope,  on  which  this  analy¬ 
sis  is  chiefly  based,  results.  Further  details  on  the  Gaussian  spectrum  are 
given  in  Appendix  A. 

B.  MODEL 

Since  a  one -dimensional  model  of  the  electron  beam  is  used,  a  uni¬ 
velocity  electron  beam  is  incident  at  the  entrance  plane  z  =  0  ,  moving  in 
the  +z  direction  in  confined  flow.  This  assumption  of  a  very  strong  longi¬ 
tudinal  magnetic  field  will  depress  the  potential  at  the  center  of  the  beam 
so  that  peripheral  electrons  travel  faster  than  axial  ones,  introducing  a 
phase  difference  between  the  radio-frequency  current  carried  by  different 
beam  segments.  This  difficulty  is  overcome  by  assuming  the  existence  of 
a  thread  of  positive  ions  along  the  axis  of  the  beam,  just  sufficient  to  neu¬ 
tralize  the  charge  density  of  electrons;  thus,  the  effects  of  depressing  the 
potential  across  the  beam  caused  by  space  charge  are  neglected,  and  so 
also  are  variations  in  electron  velocities  caused  by  thermal  noise.  Elec¬ 
tron  velocities  are  assumed  small  compared  to  the  velocity  of  light,  per¬ 
mitting  a  nonrelativistic  treatment  of  the  problem.  In  the  analysis  of  a 
double-frequency  signal,  however,  the  change  in  response  caused  by  rela¬ 
tivistic  effects  is  studied.  The  electric  field  is  assumed  constant  through¬ 
out  the  cross  section  of  the  klystron  beam. 
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II.  RESPONSE  OF  KLYSTRON  TO  DOUBLE -FREQUENCY  INPUT 


In  this  section,  the  response  of  a  klystron  to 

input  will  be  treated.  The  signal  is  V  =  Vj  sinto^t 

^o!  1^1 

a  double -frequency 

+  sinoj^t.  The 

Azl  1^3 
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four  planes  of  reference  in  a  two-cavity  klystron  are  represented  by 
,  A^  ,  Aj  ,  and  the  subscripts  o,  1,2,  3,  respectively  will  be 
used  to  identify  quantities  in  the  respective  planes. 

A.  FIRST  GAP 

Applying  Newton's  second  law  of  motion,  we  have 

,2  eV,  eV, 

-  =  -  sinw.t  +  -  sinw^t  ; 

dt2  md^  i  mdj^ 

therefore 

-eV,  eV 

—  =  -  cosw.t  -  -  cosupt  +  k  , 

dt  md^Uj  ^  1 
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where  k,  is  a  constant  to  be  determined  from  the  initial  condition:  at  t  =  t  . 

1  o 

the  velocity  dz/dt  =  ;  therefore 


dz  _ 

eV^ 

M  + 

cos  w,  t 

1  o 

cos  w^t 

COS  OJot 

Z  0 

cos 

dt 

°  md^ 

.  “l 

U) 

1  J 

mdj 

L  2 

“2 

(2.  1) 


Integrating  again  gives 


u  t  + 
o 


e  V 


md 


cos<(j,t  eV  cosuj_,t 
1  o  ^  2  2  o 


1 

e  V, 


<ji 


md 


1 

e  V. 


- —  sinw.t  + 

md^cj^ 


md,  w 


sinoj^t 


I'^Z 


+  k_ 


where  is  a  constant,  evaluated  by  inserting  the  initial  condition:  at 
t  =  t^  ,  the  distance  z  =  0  ;  therefore 


eV,  cosw.t  e  V, 

,  1  1  o  ,  2 

u  +  -  -  +  - 

°  md,  w,  md,w 


1 

e  V 


1 


cos  co-,t 
2  o 


e  V, 


(t  -  t 


—  (sinw,t  -  sinw.t  )  -  -  .(sinw-,t  -  sin.w,t  ) 

2  1  1  o'  ,2  2  2  o' 

md,w^ 


Putting 


Vl  _  ^2  _  _ 

V  ’  V  °'Z’  U  '^1  ’ 

O  O  O 


=  'f'l 


gives 


z  = 


a  u  a,u 

,  1  o  ^  ,  2  o 

u  +  -  cosw.t  +  -  cosu„t 

o  Tii  io  T.'i  2o 


Z<^ 


1 


2(j> 


1 


(t  -  t  ) 
'  o' 
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2cj> 


u 

Z  o 


(sinw,t  -sinw,t)  4-  -  (sinw  t  -  sinca^t)  .  (2,2) 

iO  1  oi**  Z  O  L 

2(j>i«2 


The  distance  is  z-d  at  t  —  t^  ,  where 


h  "  *0  ’ 


di  6  ^ 

r,  =  transit  time  in  first  gap  =  —  +  — 

1  ®  ^  u  u> 

o  o 


5  =  correction  factor, 


Wi  +  ^2 

01  =  - 

°  m  +  n 


o>-),  =  no)  . 
2  o 


m  and  n  being  two  numbers  defined  by  -w^  =  mo>^  and 


Substitution  of  this  condition  in  Equation  (2.2)  gives, 


/  ""l  “^2  \ 

d  =  u  1  + — -  cosw.t  +  -  cos«,t  r, 

°  V  °  2*1  7 


°'l'^o  '^2'^o 

+  -  (sinoj.t  -sinco.t)  +  - rr -  (sinw^t  -  sino).,t) 

-,,1  '  lo  1'  2o  2' 

2c|)^oj^  2cj>jW2 


•  (2.3) 


The  following  assumptions  are  made; 


i  j.  '  1  c  )  sin  It.,' 

o>,t  -cb,  -  —  0,  =  Uo,t  -  9, 

1  ^1  «  1  )  cos  \  1  ^1 


(0),  ' 

j.  j,"  2  - 

co,t  -  cb  -  —  6, 
2  ^1  0)  1 
o  , 


O  III  .  ,  »  ' 

o>,t  -  4), 
cos  '  2  ^1 
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The  products  <^2^1  negligible,  and  Equation  (2.3)  reduces  to, 


u  5 , 
o  1 

u 


ad,  a  ^  d ,  a  u 

- r-  cos  («,t  -tf>  )  +, — — cos  (wpt  -  cj)'')  +  — 

2=1^2  ■  ^^1^2 


^2^  n 

+  - fsinlw^t  -  .  sinw2tl 

7.rf>  L  J 


j^sin  (Wjt  -<t>  2)  -  sincojtj 


2c|>^u>2 


°’l'^l  I  r  “Z'^l 

=  - —  (cos  w,t  cos  (jil.  +  sinw.t  sinct>  )  + - -  (cosw^t  coscf)  +  sinw,t  sin<j)'') 

■  T  ,  '  '  1  .  1  1  1  T  j, "  ^  -l  ^  -i 


Zcj), 


.  + 


2tj> 


1 


a ,  u 
lor. 

_  SI 


a^u 
2  o 


sinoj^t  cos  <|>2  -  cos  u>2 1  sin  4>|  -  sinWjtj 
sinw^t  cos  -  costo^t  sin  -  sinui^tj 


therefore' 


6,  = 


a. 


2  ■■ 

2ct)''n 


1^(1  -  cos  <|)|  -  <j>^  sin  <^|)  sinW2t  +  (sin  <132  -  ^  cos  (jjj)  cosWjtj 


—  i  1  **  •  "  I  /  •  1"  x'*  x"\ 

+  - PT^  1(1  -cos  9^  "^1  sin<P'^)  sinco^t  +  (sin  9^  -  9^  cos  92)  cosw^tj 

(2.4) 


Thus,  the  expression  for  the  correction  factor  as  given. by  Equation  (2,4) 
is  a  superposition  for  each  frequency,  considered  separately.  The  induced 
current  is  calculated  by  Ramo's  theorem;  i.  e. 


e  Vt  /coscoTt  coswTtN 

+  -4  ( —5^  -  -^n  di 


mdi  I 


u  r,  + 


1  u 

1  o 


r  '  2^: 


sin^Oj^t  -  sinajj{t  -  F^) 


0), 


-  cosw^t 


2^0 

24;'  L 


sincj^t  -  sinw^Ct  " 


w- 


cosw^t 


After  simplifying,  and  in  the  process  neglecting  m5 j  cosw^^  cosWj,t, 

'  we  have 


i,  =  I  +  I  a, 
1  o  o  1 


^1  -  cos  4>2  sintjj^X  /sin<)>2  1  +  cos  (jj^ 

sinw, t  +  I  -  -  - : - 1  cos w , t 


2  c|) 


1 


2cj> 


'1' 


1 


■  2 
1 


2(j)j 


'1‘ 


+  I  a, 
o  2 


/ 1  -  cos  (J)^  sin  (j)^ 


2<)) 


Z<|) 


sin  cj)j  1  +  cos  (f>j 

..  — 2 - r:? — -} 


<t> 


1 


2({> 


(2.5) 


'From  Equation  (2.  5),  it  is  evident  that  the  effects  of  the  different  frequency 
components  of  excitation  on  the  induced  current  are  mutually  independent. 
The  double -frequency  case  can  therefore  be  extended  to  the  multiple-fre¬ 
quency  input,  and  it  can  be  concluded  that  the  different  frequency  effects 
are  independent  of  each  other,  subject  to  the  approximations  made  in  this 
section. 


B.  DRIFT  REGION 

Thus  the  electrons  enter  the  drift  space  with  both  velocity  and  current 
modulation.  From  Equation  (2.  1), 
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e 

u,  =  u,  + - - 

^  °  md,a) 


1  I 
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I  1  U  0)  11 

\  o  o/ 


_!Zi_ 

mdico^ 


f 

cos  COo  t, - - 

1  1  U  OJ 
\  o  Oj 


-  cos 


Resorting  to  approximation  made  in  Section  A  gives 


sin  <|>2 

u,  =  u  +  —  u  a, - - — 

1  o  ^  o  1 


sin  cf) 


1  \  1  ? 

sin  jw,t,  +  —  u  a,  - 

111  Z  Z  o  Z 


‘t’l 


sin  WTt, - 

■  2  1  2 


Putting 


(2.6) 


sin  A, 


sin  4) 


where  and  are  defined  as  the  gap-coupling  coefficients  for  the  re¬ 
spective  frequency  components  of  excitation,  gives  the  tim6  of  arrival  at 
plane  as 

S, 


t^  -  t^  + 


1  +  j  o-iPi  smlwiti-— j  +j°-2^1  ®^”'(‘^2^l‘  “ 


For  small  excitations,  we  get 


t,  =  t,  +  — ; 
2  1  u 


a^pj  sinLt^  --1)  -ia^p;'  sin  Ltl  '  — 


I 
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I 


c^iaiPi 


COS 


o  ' ' 


COS 


The  equation  for  conservation  of  charge  gives, 


.  _  h 

1 

I 

0 

^2 

dt^ 

dt^ 

dt^ 

1 

1 

|t3 

dtj^ 

dt 

0 

Again  t,  =  t  +r,=t  + 
®  1  0  1  0 

li 

u 

0 

.!i 

CO 

0 

;  therefore 

dt. 


dt^ 


1+-^  fill 

dto 


as  is  independent  of  ,  and 


dt  a 

_L  =  1  +  —L 

2cj, 


24,, 


-  <{)j  sin  cl)|  -  cos  <}>^j  cos  +  (j>jj  -  ^sincj)^-  4i|  coscji^^j  'sin^Wj^t^+  cj>^j 


^0)^*0+ ‘t’l)  -  (sincj)^'-  cj)^'  coscjjj  jsin 


sincjjj  -  cos4>j  j  cos 


As  a  result. 


i,  =  I 
2  o 


s!a,p: 


1  cos  Kt^  -- 


f ,  h'' 

-  cos  (Wot,  -  - 

2  \  2  1  2 


^  +  _i_  |i^l  -  cj>j  sin(|)^  -  cos  (f)^j  cos |w^t^  +  4> -  ^sintj)^^  -  (j)^  cos<t>^j'  Sin^w^t^+4)' 
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2^,,  f  "'*>1  sin'll  -  cos<j);'j  cos  -  (sin<j)”-c|,J'  coscj>”j  sinjo^t^  +  cj,;' 

(2,8) 


where 


With 


S  = 
1 


“l^l 


w.,S , 

S"  =  ^  ^ 


U 


o 


sin 

CO 


s  (“I'o  -  ♦;)  =  M 


and  with  approximations  similar  to  those  made  previously,  the  current  i^ 
can  be  expressed  as  a  function  of  t^  ,  and  hence  t^  , 

The  periodicity  of  i^  is  obvious,  and  therefore  i^  can  be  expanded 

1  _ -y  ^  j.  ^ 

2 


in  Fourier  series  as  follows.  First,  t^  will  be  related  to  t^  ,  as  follows 


d  S  . 

t,  =  t  +  -i  + 

2  o  u  u 

o  o 


1.  e.  ,  the  term  is  neglected,  while  the  approximations  with  regard 

to  the  sinusoidal  terms  are  justified:;  then 


It  —  *—  a  -|- 
2  2  0 


r--l 


cosn^wt^  -  -  (j)J  ^ '^r  ^^osn^wt^  -  -  (fij 


r=l 


where 


to  S, 
o  1 


4)1  =_ 


a  =  21 
o  o 
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and 


^  T 

=  -  f  — 

dt^ 


TT 

=  -  f  1 

ir  /  o 

-IT 


COS  r(a)t2  "  "  ‘I’j)  ^(wt^) 


cos  r 


oj  t  -k,  sin  mw  t  +  — —  i  -  k,  sin  nw  t  +  — 
ool  100  2/1  oo 


5 

2 


1  fi  .  r ,  ,•  .  /  , , u'  .  f 

=  —  /I  smrojt  -  k,  sin  mco  t  +  —  -  k.  sin  nco  t  +  — 

TT  J  O  OO  i  \  O  O  I  1  l002 


-ir 


d(aj  t  ), 
o  o 


d(co  t  ), 
'  o  o'  ’ 


(2.9) 


where 


k,  = 


SiaiPi 


,n 

"^1  =  -2  — 


The  coefficients  and  are  simplified  by  a  method  indicated  in  Appen- • 
dix  B. 

The  general  equation  of  motion  has  to  be  formulated  for  the  calcula¬ 
tion  of  the  induced  current  in  plane  A^,  and  the  method  adopted  is  similar 
to  that  used  in  the  first  gap,  with  only  the  initial  conditions  different.  This 
is  done  later  for  the  more  complicated  Gaussian  spectrum  and  will  be 
omitted  here.  The  main  purpose  of  this  section  is  to  show  the  absence  of 
intermodulation  of  the  different  frequency  components  at  the  output,  subject, 
of  course,  to  the  approximations  made.  It  must  be  stated  that  the  preceding 
analysis  was  based  upon  frequencies  and  not  being  very  far  from 

the  central  frequency  in  the  frequency  spectrum,  i.  e.  ,  the  numbers  m 
and  n  should  not  be  much  greater  than  1.  As  the  main  purpose  of  this  study 
is  the  extension  of  this  analysis  to  the  response  of  a  klystron  to  short  pulses 


-12  - 


with  a  narrow  bandwidth  of  frequencies,  the  assumptions  made  are  com¬ 
patible  with  the  condition  desired. 

No  account  has  been  taken  of  the  relativistic  variation  of  mass  with 
velocities.  This  problem  becomes  especially  serious  when  the  beam  voltage 
is  large  in  high-power  klystrons  and  where  the  very  hard  X-rays  produced 
present  an  additional  hazard  to  the  operating  personnel.  A  simple  treatment 
will  be  given  of  the  relativistic  effects  on  the  response,  using  the  same 
■model  as  before. 

According  to  Einstein,  nothing  can  move  with  a  speed  greater  than 
the  speed  of  light.  Newtonian  mechanics  comisined  with  this  postulate  de¬ 
mands  that  a  mass  subjected  to  a  constant  force  must  be  accelerated  till 
the  speed  of' light  is  attained;  but,  as  the  force  is  still  present,  the  speed 
must  still  increase,  which  is  impossible.  This  ambiguity  is  solved  by  ac¬ 
cepting  the  increase  of  mass  with  velocity,  and  assuming  that  mass  is  a 

manifestation  of  energy,  the  two  related  to  each  other  by  the  famous  equa- 
2 

tion  w  =  c  m  ,  where  c  is  the  velocity  of  light. 

2 

An  increa-se  in  mass,  dm,  when  accelerated, results  in  c  dm  = 

=  dw  =  F  ds  ,  where  F  is  the  applied  force  over  a  distance  ds  .  Newton's, 
second  law  gives 

F  =  ^  fmv)  ; 

therefore 

=  / 

Equating  the  integrands  and  separating  variables,  we  have 

dm  V  dv 


—  (mv)  ds 
dt 


r 


vd  (mv) 
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Assuming  that  rest  mass  equals  m^  .  gives,  by  integration, 

m 

°  (2.10) 

the  equation  that  demonstrates  the  variation  of  mass  with  velocity.  In  this 
case,  Newton's  force  equation  gives 


With  u(z,  t)  =  +  v(z)  e'^^^  ,  where 

u(z,  t)  =  total  electron  velocity, 

u^(z)  =  d-c  beam  velocity, 

v(z)  =  amplitude  of  a-c  velocity. 


The  basic  assumption  will  be  that  v  «  c,  which  is  justified,  since 
the  signal  voltage  is  not  sufficiently  high  in  practice  to  make  the  a-c  velo¬ 
city  appreciable  in  comparison  to  the  speed  of  light.  Using  the  Taylor 
series  expansion,  we  have 


14:-. 


therefore 


_d_ 

dt 


u 


ve 


jwt 


2\1/2 


m 


(2.  11) 


u 

II  -  — 

c 


As  the  excitation  is  variational,  separating  Equation  (2.  11)  into  the  d-c 
and  a-c  parts  gives, 

=  0 


sinw^t  +  sinoi^tj  ; 


(2.  12a) 


therefore 


V 


1-a-c 


u 

[1  -  -H- 


3/2 


m  d, 
o  1 


cosw^t  cosw^t 


u 


1 


w- 


+  k. 


where  k_  is  a  constant  evaluated  from  the  initial  condition  that  at  t  =  t  , 
3  o 

V ,  =  0  .  Finally 

i  3.  —  C 


1  a-c 


/  2i\ 

r 

fi-V 

1  e 

V  1 

( cos  w,  t 
'  1  o 

cos  w^t  \ 

V 

1  m  d 

o  1 

^  '^1 

“1  / 

+  V, 


/  cos  w-t 
/  2  o 

cos 

\  “2 

•  “2  j 
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The  total  velocity  at  any  instant  in  the  gap  is 


dz 

=  u  +  V, 

dt  o  1  a-c 


Integrating  again  and  using  the  initial  condition  that  z  =  0  at  t  =  t^  , 


we  have 


z  = 


3/2 


u 


u  +  1  -  _ 

O  I  m  d 

*=  /  o  1 


/V,  cosWit  V,  cosoj^t 

e  /  1  L2.  +  .  ^  ^  ° 


a> 


1 


iJi- 


(t  - 


/  2V/‘- 

/  Cl  \ 

/ 

11  -  — 

1  ®  ( 

1  _2  , 

m  d  \ 

V  ^  y 

o  1  \ 

V,  sinw,t  V,  sinco.t  V-,  sinw^t  sina>-,t  ^ 

1  l_o  _  1  1  o  ^  2  2  o  _  _2 _ 2  o 


cot 


Now  z  r=  d  at  t  =  tj^  =  t^  +  .  Following  the  procedure  of  Section  A  of 

this  chapter,  and  making  similar  approximations,  we  obtain  the  expression 
for  the  correction  factor  in  this  equation: 


«1  = 


^1^1 

2ct)^'w^ 


1  -  cos  sin<j>j^'j  sincojt  +  ^sin  cos  cos  w^^t 


2<j) 


CO  r 

o  I  ,  XM  i”  .  X »»  .  ^  ,  II  .  t» 

- -  1  -cos9^  -  <P^  sin<}>j  sinw^t  +  (sincj)^  -  c)>j^  cosc|)^ 

I'w-,  L  ' 

ZnV2 


I'Z 


cos  W  jt 


u 

■o 


(2,  12b) 


Equation  (2.  12b)  is  very  similar  to  Equation  (2.  4)  and  shows  that  the 

/  7n3/2/2 

correction  factor  is  only  multiplied  by  a  constant  (1  -Uq)  /c  when  rela¬ 
tivistic  effects  are  included.  Again,  the  total  current  induced  as  a  result  of 
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the  passage  of  electrons  in  the  interval  t^  =  t  -  ,  and  =  t  is 

I  * 


i  -  -2^  f  u  dt 

1  d  J  o 

t-r. 


1 

Carrying  out  this  integration  as  in  the  nonrelativistic  case,  we  have 


sinw^t  + 


sintjjj^  (1  +  coscj)^) 


2c)>. 


2(j) 


|cos  Wjt) 


tcos4>'  sin<j)''\  /sincj)''  l+costjj'' 

1 - ^ - i  sinw,t  +  ..—  1 - 1. 


cosw^t 


(2.  13) 

From  Equation  (2.  13)  the  relativistic  effects  on  the  induced  current 

in  the  first  gap  are  very  clearly  observed..  The  explicit  effect  on  the  r-f 

current  resulting  from  electrons  at  high  beam  voltages  subject  to  the  simple 

approximations  made  should  be  noted.  For  a  particular  beam  voltage,  the 

2\  3/2 


r-f  induced  current  is  lower  by  the  factor 


fl-u^ 


when  the  rela- 


I 

tivistic  variation  of  mass  with  velocity  is  taken  into  account.  A  curve  has 
been  plotted  to  show  this  effect  (Figure  2), 


For  low  voltages  and  hence  low  values  of  u^  ,  the  factor  ^1  -  j 
is  equal  to  unity,  and  the  result  becomes  similar  to  that  derived  for  the 
nonrelativistic  case.  For  u^/c  =  > 


3/2 


vV2 


7  - 


=  0.999 


'4 
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10 


Figure  Z 


BEAM  VOLTAGE  EQUIVALENT  OF  VELOCITY 


10 


6 


Graph  Showing  Relativistic  Effects  on  Current  Response. 
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III.  BALLISTIC  ANALYSIS 


A.  FIRST-ORDER  BUNCHING  THEORY 

The  problem  under  consideration  is  formulated  as  follows:  A  beam 
of  parallel  electrons  which  have  been  accelerated  through  a  potential  of 
volts  is  passed  through  the  grids  of  a  resonator  across  which  there  appears 
a  voltage  sinwty,'.  The  resultant  electric  field  is  assumed  paral¬ 

lel  to  the  electron  motion.  Since  the  velocity  of  an  electron  is  proportional 
to  the  square  root  of  the  voltage  through  which  it  has  been  accelerated,  the 
velocity  with  which  an  electron  emerges  from  the  firsti/,  or  bunching,,  resona¬ 
tor  of  a  two -resonator  klystron  will  be^ 


V  = 
a 


'1  + 


Pl^l 

V 


sinwt  , 


■=  V 


1  + 


sinwt 


I  2e 

where  u^  =  /—  is  the  d-c  beam  velocity,  is  the  gap-coupling,  co¬ 
efficient  (taking  into  account  the  effect  of  the  gap  transit  angle),  and  = 
V^/V^  .  Here  may  not  be  related  to  the  gap  transit  angle  in  the  same 
way  a§  for  the  sinusoidal  case,  but  it  obeys  the  general  definition  of  the 
ratio  of  the  velocity  gained  in  the  real  gap  with  across  it  to  the  velocity 

gained  in  an  infinitely  narrow  gap  with  across  it  and-,  as  such,  is  always 

less  than  1.  The  time  taken  by  an  electron  to  move  a  certain  distance  along 
thp  beam  depends  upon  the  point  on  the  cycle  at  which  it  passed  through  the 
resonator  gap  as  well  as  upon  the  magnitude  of  the  gap  voltage.  If  is 
the  drift  length,  t^,  the  tipae  at  which  the  electron  leaves  the  first  resona¬ 
tor,  and  t^  the  time  of  arrival  at  the  catcher;  then 
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t  + 
o 


-at 


TJI 


I  1  +  e 


sinwt 


Now,  if  the  modulation  factor  is  small  compared  to  unity;  then  the  fol¬ 
lowing  approximation  is  reasonably  valid: 


t 


2 


t  + 

O' 


u 

o 


To  find  the  current  associated  with  the  electron  bunches, one  must 
remember  that  the  principle  of  conservation  of  charge  applies  to  electron 
bunches  for  an  interval  with  corresponding  departure  and  arrival  times.  The 
electron  stream  is  subject  to  the  conservation  of  charge,  so  that 


I  dt 
o  o 


where  i2  is  the  catcher  current;  hence 
I 


^2  ■ 


dt^ 


dt 


1  + 


u  2 

O  ^ 


- 

/  2  2  \ 

f  -at^  -at^  ^ 

-  2  at  e  sin  wt  +  a)  e  coscDt 

O  0  o. 

V  / 

-  a9^/w^  a9  -a.Q^ /<j)^ 

1  -  k  ( e  °  cos  9  -  2  — 2.  e  °  sin  9 

o, 
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As  this  expression 


where  the  bunching  parameter,  —  — -  =  S,  - — 

is  aperiodic,  it  cannot  be  represented  by  Fourier  series,  contrary  to  the 

2 

analysis  carried  out  by  Beck  for  the  pure  sine  wave.  Since  the  catcher  re¬ 
sponse  is  desired,  it  has  been  found  convenient  to  plot  the  output  current 
versus  t^  for  different  bunching  parameters,  similar  to  the  treatment 
given  by  Spang enburg^  for  the  sinusoidal  excitation.  Since  i,  —  f(0  )  = 

=  f(-9  )  the  curves  are  expected  to  be  symmetrical.  If  the  Gaussian  spec- 

°  -at^ 

trum  is  represented  by  Ve  cos  wt,  then 


^2  " 


1  +  k 


°  sin  0  +  2 — 2.  e  °  cos  0 

°  0,2 


In  this  case  i^  =  f(0  )  ^  f(-0  •) .  For  negative  values  of  0  ,  we  have 


^2  “ 


o 


1  -  k 


-a0  /o,  a0  a0  ,/o, 

o'  •  Q  ,  o  o  o' 

:  sin  0  +2  — —  e 

o  2 

CO 


cos  0 


The  choice  of  a/co  is  governed  by  the  following  consideration.  The 

-at^ 

envelope  has  its  maximum  value  at  t  =  0  (since  the  envelope  is  V  e  cos 
cot)  and  is  supposed  to  fall  to  l/e  of  its  maximum  when  cot  =  lOir,  so  that 
10  r-f  cycles  are  enclosed  between  the  points  where  the  amplitude  is  l/e 
of  the  maximum  value;  therefore 


4-  (10^)^  •=  1 


a 

T 


CO 


(1017)“^ 
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Graphs  have  been  drawn  for  the  output  current  versus  exit  tinae  for  differ- 

o  ■j-Z 

ent  values  of  bunching  parameter  for  the  envelopes  V^e"  ''  sinwt,  and 
2 

e"^^  coscot,.  In  the  first  case  the  infinite  peaks  occur  for  values  of 
9^  satisfying  the  transcendental  equation, 


The  values  of  6^  at  which  infinite  peaks  occur  are  found  graphically;  For 
k  <  1  ,  there  are  no  infinite  peaks,  as  is  evident  from  the  equations  also 
(Figure  3);  for  k  =  1  ,  one  infinite  peak  occurs  (Figure  4);  for  k  =  1.5  , 
there  are  14  infinite  peaks  (Figure  5);  for  k  =  2  ,  there  are  18  infinite  peaks. 

This  can  be  justified  as  follows.  In  the  pure  sine -wave  case,  two 
infinite  peaks  occur  for  k  >  1  ;  hence,  for  simplicity,  we  associate  two  in¬ 
finite  peaks  with  two  peaks  of  the  excitation  signal.  In  the  Gaussian  en¬ 
velope,  the  seventh  peak  occurs  on  either  side  of  t  =0.,  when  wt  =+  6.5tt. 

-at^ 

When  cjt^  =  lOu,  then  e  =  l/e;  therefore  when  cot^  =  6.5it  ,  then 


1 


0.425 

e 


The  response  for  k  =  2  has  not  been  actually  plotted,  for  it  is  not  ex¬ 
pected  to  be  very  dissimilar  from  the  k  =  1.5  response  plot.  The  num¬ 
ber  of  infinite  peaks,  however,  were  determined  by  finding  the  points  of 
intersection  of  the  curves: 


a9 


a9 


cos  9  -  2 

o 


sin  9 


and  —  e 
2  ■ 
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Figure  4.  Output  Current  versus  Exit  Time  wt^  -S  . 
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Figure  s,  Output  Current  versus  Exit  Time  wt^  -  S  . 


Figure  6,  Output  Current  i^  versus  Exit  Time  wt^  -S. 
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The  value  of  k  when  u)t  =  6,5it  (if  k  =1.5,  when  t  =  0)  is 

o  o  o 


1.5 


6,5 


,0.425 


1.5 

1.5 


=  1 


Hence,  14  peaks  are.  enclosed  in  the  region  for  which  k  >  1  . 

Again,  the  ninth  peak  occurs  on  either  side  of  t^  =  0 ,  when  wt^  =  8.5'iT, 

so  that  this  region  encloses  18  peaks.  The  value  of  k  at  wt^  =  8,  Sir  (if 

k  ■=  2  ,  at  t  =  0)  is 
o  ’  o  ' 


1 


8-5  "  gO.725  2.02  — 


Thue,  the  occurrence  of  18  infinite  peaks  for  k  =  2  is  justified. 

-at^ 

If  the  envelope  is  e  cos  ut ,  no  infinite  peaks  occur  for 

k  =  1  (Figures  6- and  7').  For  k=  1.5,  there  are  12  infinite  peaks  (Figure8  ), 

-at^ 

This  has  also  been  justified  by  a  process  similar  to  the  V^e  sin  (ot 
case. 


B.  KLYS.T,RON  RESPONSE  WITH  VARIABLE  GAP  LENGTH 
From  Newton's  equation  of  motion. 


A 

dt 


e  V 


1 


-at  , 

e  cos  wt 


hence 


^  =  f  ^ 

dt  j  md 


-at^ 

e  cos  a>t  dt  +  C 


where  C  is  a  constant,' When  t  =  t  , '  then  dz/dt  =  U  .  Thus, 

,1  o  '  o 
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Figure  8. 
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IOtt 


-26- 


cos  wt  dt  =  C 


t 


from  which 


<i-  .  r  ■“'"i 

—  =  u  +  /  -  e  cos  wt  dt  ; 

dt  °  j 


hence  the  exit  velocity  from  the  first  gap  at  t  =  is 


u,  =  u  +  / 
1  o  J 


e  V  Z 

1  -at  .  ,, 

-  e  cos  wt  dt 

md 


The  time  of  arrival  at  the  catcher,  if  S^.  is  the  drift  distance,  is 


t,  =  t  + 
2  o 


r‘  -Z 

u  +  -  /  e  CO'S  wt  dt 


o  md 


therefore 


Si  /  eV 

.  t,  ■=  t  +  — L  1  + - — 

2  o  ^ 


e  COS  wt  dt 


o 


Here  an  approximation  will  be  made  to  permit  analytical  computation: 
t^  =  t^  +  r  ,  where  the  transit  angle  r  is  (^/^q)  +  (^/w)  • 

We  assume  here  that  6  =  0,  so  that  r  =  d/u^  .  Again  for  small 
modulation,  the  following  expansion  is  permissible,  at  least  for  a  first  ap¬ 
proximation: 
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t 


I 


\ 


dt^  S  ^  ® 

dt  ^  ^  md 


-a(t^+r)^  -at^ 

-e  cosw(t  +  r)  +  e  cos  wt 

o'  o 


Si  ^ 

md 


‘I'o 


wh 


ere  (j>^=wd/u^,  and  k'^Spi/Z,  hence, 


^  =  1  t 


o  L 


_  a  .2 

”  2  o 

cos  (9  +(b  )  -  e  ^  cos  0 
o  ^o  o 


thus 


^2  “ 


dt2 

dt„ 


1  - 


-  (9  +cj) 

2  '  o  ^o' 

e  cos  (0  +  *  -e 

o  ^o 


a  -2 
2  °o 


cos  0 


oj 


The  above  expression  for  ,  though  approximate,  does  give  an  idea  of  the 
catcher  response  at  least  fpr  small  signals i 


The  expression  for  i^  becomes  a  poor  approximation  for  large  gap  angles 
also.  The  graph  of  i^  versus  0  with  tj>  =  tt/Z  has  been  plotted  (Figure 
10),  taking  k  =  0.5,  and  1.0;  and  tSe  differe^ice  between  this  and  the  first- 
order  bunching  theory  is  obvious.  It  has  been  found  convenient  to  plot  i^ 
versus  0^  ,  in  this  case,  rather  than  i-  versus  0^  ,  as  we  are  only  in¬ 
terested  in  an  approximate  estimate  of  the  response  curve. 
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Figure  9.  Output  Current  versus  Exit  Time  wt^  -  S  . 


'  i 


Figure  10.  Output  Current  i,  versus  0 


C.  RESPONSE  WITH  FIXED  ENVELOPE  AND  VARYING  FREQUENCY 

Using  the  first-order  bunching  theory,  and  taking  k  =  0.5  ,  we  have 
plotted  the  catcher  current  response  for  a  fixed  envelope  for  three  differ¬ 
ent  r-f  frequencies  (Figure  9)-  It  has  been  found  that  the  shape  of  the 
response  envelope  is  almost  independent  of  the  change  in  frequency,  es¬ 
pecially  at  high  values  of  . 
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IV.  LARGE  SIGNALS  AND  FINITE  GAPS 


I . 

'  The  ballistic  analysis  for  large  signals  with  finite  gaps  is  carried 

!  out  as  follows..  Analysis  is  based  on  a  two-cavity  klystron  (Figure  11). 
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Figure  11.  Schematic  of  Model  for  Velocity -modulated  Tube. 

In  the  diagram,  A^,  A^,  A^,  A^  represent  the  four  planes  under  con¬ 
sideration.  The  numbers  0,  1,  2,  3  will  be  used  to  identify  quantities  in  the 
respective  planes. 


A.  FIRST  GAP  REGION 
1.  Induced  R-F  Current 

Confining  our  attention  first  to  the  motion  of  an  electron  in  the  first 

2 

gap,  excited  by  a  Gaussian  pulse  V^  e"^'"  coswt,  and  applying  Newton's 
second  law  of  motion,  we  get 


d^z 


dt 


Hi 

md 


-at 


cos  wt 


I 


-31- 


therefore 


where  K,  =  eV,/md,  and  C,  is  a  constant.  Now  at  t  =  t  ,  dz/dt  =  u  , 
1  r  1  1  o  '  o  ’ 

from  which  is  evaluated;  thus 
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Integrating  again,  we  get 
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where  =  constant.  This  constant  is  evaluated  by  using  the  condition 

that  at  z  =  0  ,  t  =  t  ;  therefore 
o 
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o'  o'  2 
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When  t  =  tj^  ,  z  =  dj  .  Also  the  transit  time  is  =  t^  -t^  ;  therefore 
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Suppose  F  =  —  +  — ,  where  6  is  the  correction  factor;  then  wd,/u  = 
X  u  Qj  1  1  o 

o  ^ 

=  4>2  is  the  d-c  transit  angle.  From  Equation  (4.  3), 
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K-n/tT  -«  /4a 
1  e 


2\/a 


(4>2  +  Sj) 
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erf's/a  t  +— 7=-  /  e  °  sin'v/a’ t  y  dy 


O  ^/^]• 


Writing  as  t  and  replacing  by  tj  -  ^l"u“  ~ 


In  Equation  (4.4),  5^is  implicit,  and  approximation  must  be  made  to  obtain 
an  explicit  expression  for  it. 
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We  shall  assume  that  5^  is  small  so  that 
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Now  as  is  small  compared  to  t,  the  following  approximation  is  justi¬ 
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Substituting  these  in  Equation  (4.4),  we  have 
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u 
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This  equation  obviously  gives  a  quadratic  in  6^  ,  which  might  be  complica¬ 
ted  to  solve.  For  a  simpler  solution,  we  will  neglect  the  5^  term.;  there- 
fore 
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(4.5) 


Thus  Equation  (4.  5)  gives  an  explicit  expression  for 
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According  to  Ramo's  theorem,  the  current  at  time  t  resulting  from  the 


charge  entering  the  gap  in  the  interval  dt  between  time  t  and  t,  +  dt  is 
&  to&r-  Q  OOO 


I 


di ,  =  -  u  dt 


where  u  =  velocity  at  time  t.  The  total  current  induced  is 


i. 
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Substituting  Equation  (4.  lb)  in  the  preceding  integral,  we  have 
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2  Ft 


'  /  *2  Z> 
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-  e  erfi^Jk  t  dt  + -  /  e  °  sin'\la  dy  dt^ 


(4.7) 


The  above  expression  is  considerably  simplified  if  the  following  approxi¬ 
mations  are  accepted  as  valid; 
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This  is  very  nearly  true,  especially  for  small  gaps;  then,  Equation  (4.  7) 


reduces  to 
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It  is  interesting  to  note  in  this  case  that  the  time -dependent  component  of 
current  i^  is  directly  proportional  to  the  correction  factor  5  , 

A  more  accurate  simplification  of  Equation  (4.  7)  than  that  given  by 
Equation  (4.8)  is  obtained  if 
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Putting  this  in  Equation  (4.  7),  we  have 


i,  =  I  +  — 

1  o  d, 


u  5 , 
o  1 


K, 


-ifi. 

2^y^| 


After  expansion  in  a  Taylor  series,  we  get 
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If  the  5^  term  is  neglected, 
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m  2 

where  V  =  d-c  beam  voltage  =  —  u 
o  2e  o 


2.  Velocity  Modulation 

From  Equation  (4.  la),  we  have 
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Making  suitable  approximations  to  simplify  solution  gives 


^  2n^ 


2 

nAt 


e-‘"  dt 


^Ar  I  % 


-  f^ia  t,  - 


^^t  -iii- 
°  2'\^ 


from  which 


V 


K, 


_J_  ^  nAt 

u  ^  2  u 

o  o 


'di  6^Y  -at^ 

—  +  —  e 
u  to  / 

O 


cos  tot. 


1 


(i.  10) 


This  is  justified  if  the  gaps  are  very  short  and  the  transit  time  small.  For 
larger  gaps,  a  better  approximation  would  be- 


^^t  ->L- 

°  2^y^ 

then 


e 


(4.  11) 


The  neglect  of  5^  in  the  exponential  is  reasonable  when  the  resulting  simpli¬ 
fication  in  computation  is  taken  into  consideration.  Note  that  when  6^=0, 
Equation  (4.  10)  reduces  to 
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1  1  1 

■  —  =  1  + -  e  cos  «t,  , 

u  2V  1 

O  o 

an  expression  similar  to  the  velocity  modulation  in  the  sinusoidal  excita¬ 
tion  case. 


B.  DRIFT  SPACE  REGION 

The  electron  beam  enters  the  drift  space  with  both  velocity  and  cur¬ 
rent  modulation,  as  shown  above,  and  it  drifts  in  a  field-free  space  resulting 
in  further  increase  in  the  harmonic  content  of  the  beam  current. 

Neglecting  space-charge  debunching,  we  have 
where  is  the  drift  length;  therefore 


t^  -  t^  + 


1 


%  ^  + 


K 


1  e  4a 


n 

erf  (' 

Na 

[  ' 

jw 


/ 


/ 

Using  the  simplified  expression  in  Equation  (4.  10)  gives 
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Using  Equation  (4.  11)  results  in  a  slight  modification: 


(4.  12) 


(4.  13) 
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This  expression  for  small  modulation  becomes 


It  is  difficult  to  differentiate  this  function  because  of  the  presence  of  6  . 
It  is  better  to  keep  the  integral  in  the  expression  for  dz/dt  rather  than 
divided  into  error  functions;  thus 
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therefore 
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But  by  Equation  (A.  1),  is  related  to  ;  therefore  theoretically,  it;  is 
possible  to  express 

i^  =  f(tj)  ,  (4.15) 

where  f(tj^)  denotes  a  function  in  t^^  .  Again  as  t^  is  related  to  tj^,  as 
shown  above,  it  is  possible  to  express  i^  as  a  function  of  t^  ,  i.  e.  , 

"2  •=  ^o  §^*2^  • 

Although  Equation  (3.  l6)  is  expected  to  be  very  complicated,  and  rigid 
mathematical  analysis  seems  highly  improbable  in  practice,  a  theoretical 
formulation  is  not  ruled  out.  It  is  desirable  to  expand  i^  in  a  series  in 
such  a  way  that  the  various  frequency  components  become  distinguishable, 
but  the  obvious  aperiodicity  of  i^  rules  out  the  possibility  of  expanding  in 
a  Fourier  series,  as  in  the  sinusoidal  case  analyzed  by  Beck. 


C.  SECOND  GAP  REGION 

At  this  stage,  attention  will  be  directed  to  the  motion  of  electrons  in 
the  second  gap.  The  beam  induces  a  voltage  on  the  grids,  and  because  the 
voltage  produces  a  change  in  beam  current  and  velocity,  we  can  consider 
this  behavior  as  a  reciprocal  relationship,  so  that  knowing  the  effect  of  the 


voltage  on  the  beam  is  equivalent  to  knowing  the  effect  of  the  changing  beam 

upon  the  induced  voltage.  With  this  in  mind,  let  us  assume  the  voltage  in- 

—  b(t“l" 

duced  in  the  'second  gap  to  have  a  spectrum  given  by  the  equation  V^e 

cos  <«)(t+  p) ,  where  p  is  a  constant  introduced  to  take  into  account  the 

-at^ 

possible  change  of  phase  with  respect  to  the  original  envelope  Vj^e  coswt. 
Although  this  assumption  might  differ  from  the  actual  physical  conditions, 
it  is  acceptable  as  a  first  approximation.  The  equation  of  motion  in  the  gap 
can  therefore  be  written  as 
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Using  this  condition  to  calculate  we  have 
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Integrating  again,  we  have 
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where  =  constant.  At  t  =  t^  ,  z  =  ;  therefore 
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Now  at  z  =  d2  +  Sj^  +  dj^  ,  we  have  t  =  t^  •  Also  t^  -  t^  =  is  the  transit 
time  in  the  second  gap;  therefore 
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where  6^  is  a  correction  factor  for  second  gap..  Then,  by  a  process  similar 
to  that  used  for  the  first  gap,  and  with  similar  approximations,  we  have 
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Thus,  the  expression  for  5^  is  very  similar  to  that  for  5^^  . 

As  before,  the  total  current  induced  as  a  result  of  the  passage  of 


charge  in  the  interval  F^  is 
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t-r. 


u  dt 

o 
dtp 
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where  u  is,  as  usual,  the  velocity  at  any  instant  t  in  the  first  gap.  As 
we  have  seen  in  the  above  analysis  dt^/dt^  can  be  expressed  as  a  function 
of  t^ ;  therefore 
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This  analysis  can  be  extended  to  n  gaps^  and  n-1  drift  spaces;  under 
these  conditions 


'dt. 


dt. 


dt , 


‘•Zn-1 


'^Z(n-l)  '^^Z(n-l)  /  '^(n.Z)'^^Z(n-Z) 


dt^ 


d^dz...dn 


u  dt 


t-r 


n 


t-r 


n-1 


t-r, 


{4.Z0) 


Expression  (4.  ZO)  is  of  theoretical  interest  since  it  indicates  the  dependence 
of  the  final  current  in  the  output  gap  upon  the  excitation  imposed  upon  all 
other  gaps. 
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V.  NONLINEAR  SPACE-CHARGE  WAVE  ANALYSIS 


Here  also  attention  will  be  directed  to  the  Gaussian  excitation, 

since  the  behavior  of  an  electron  passing  through  alternate  gaps  and  drift 

regions  with  sinusoidal  excitation  in  the  gap  has  already  been  determined. 

3 

Mclsaac  has  derived  a  general  expression  for  polarization  in  the  drift 
region.  Using  his  symbols  for  the  input  gap  and  the  drift  case,  the  polari¬ 


zation  Z^(T^  ,  T^)  in  the  drift  region  is; 


Zi(T^,  TJ  = 


0(r  + T^)  sin(Tj  -  r)  dr  . 


For  the  Gaussian  case,  6(T)  =  Ae"''*'^  , 


Zi(T^,T^)  =  - 


■Y(r  +  Tj' 


cos  (r{r  +  T^)  sin(T^  -  r)  dr 


Replacing  the  sine  and  cosine  terms  by  exponentials  and  defining 


erf  X 


.  ^  f  e-‘ 

's/ir  J 


dt  , 


we  have 


Z(T,  T  )  = 


o’  "  F  if 

(o~+i)^  r 

-je  e"'^'^  /erf  ^Jy{r+T  -erf  sjy  T  I 

\  L  °  2  V7J  L  °  2  V7J/ 


+  complex  conjugate 


(5.  1) 
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For  the  input  gap  region. 


J 

Z,(T,T  )  =  -  -L  f 
1  i  o  2D  j 


0(r  +  T^,)  sin  (Tj  -  r)  dr  ; 


0 


therefore,  for  the  case  under  consideration: 


^  j 

he  e^^  ^rf 

T  -  j 

-erf 

16  i 

Iv  1  ^ 

2'V7'  . 

(0-+1) 


-je 


/erf 

-  erf 

.  j  <‘'+‘11 

1 

2^;^. 

+  complex  conjugate 


(5.  2) 


A  method  will  be  indicated  to  express  the  complex  error  function  in 
terms  of  real  integrals.  Consider  the  expression, 

a-jb 


erf(a-jb)  = 


-z 


dz 


in  the  complex  plane  with  z  =  x  +  jy  .  As  '  e  is  an  analytic  function, 
f(z)  =  0  around  a  closed  path ofintegration.  Hence  integrating  along  the 
path  shown  by  the  arrows  gives 
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Figure  12.  Gaussian  Envelope, 


from  which 


oo 


c  !  2  ^  •  -(a^-y^)  +  2jay 

erf  (a-jb)  =  — —  /  e  dx - —  j  e  '  ^  ^  dy 


N/ir 


0 


b  b 

=  erf(x)  +-^  /  e~^^  ^  sin  Zaydy  -  /  e  ^cos  Zaydy 

J  'iir  J 


(3.  3) 


As  the  polarization  has  to  be  a  real  quantity,  and  noting  that  =  T  - 
=  T  -  =  T  -  Z  +  Zj  ,  we  have  Equation  (5,  3)  written  in  the  Z  ,  T  co¬ 

ordinate  system  as  follows: 


Zjiz,  T) 


2  Re 


_A 

16 


n/^  (r+T  -  Z+Zp 


:(or-l) 


-  erf 


>^  {1  -  Z  + 


zp  -  j 


(cr-1) 


2a/^ 


(o-+r)^ 

-Tr“e‘3T 


^y7(r+T-Z+Z^)  -  j 


(o'+i) 

2a/7 
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-erf 


n/7(t-z+z^)  -j 


(5.4) 


where  Re  represents  the  real  part.  Equation  (5.  4)  is  obviously  implicit 
in  Zj^  and  hence  cannot  be  evaluated  easily. 

Assuming  to  be  sniall,  a  series  expansion  of  the  complex  func¬ 

tions  will  be  made,  and  only  the  first-order  term  in  Z^^  will  be  taken  into 
account.  A  typical  term  is  expanded  as  follows: 


erf 


N/V  (r+T-Z+Z  )  -  j  ilzil 

A  o .  r 


=  erf 


‘\/7 (r+T-Z)-j  lllil 


+  Z^  erf 


N/7(r+T-Z)  -  j  ilzil 


Zi\f^ 


where 


erf  =  e 


'v/7 (r+T-Z)  -  j  lllil 


Z^Jy" 


Substituting  in  Equation  (5..  4),  we  have 

/( 


+  ^  Z^  e 


(<r-l)^ 
je 


'erf 


N/yir+T-Z)  -  ’ 


'^Fi(T+T-Z)  -  j 


Z\f^ 


+  \fy 


^y7  (T-z)-j 


l2 


2n/y' 


erf 


(o-+l)‘ 


Z\Pf 


-je  ^  e-J^’l 


erf 


z\F^  J 


\ly~ (r+T  -  Z) 
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-J 


.  (cr+1) 


-  J 


;  (o'+i) 


l>s!y 


+  n/y 


+  Z^e 


'^y7(r+T-Z).j  lllil 


,2^ 


2n/7 


in 


n/7  (T-Z) 


n/7 (T-Z)-j 

2a/V 


(5.5) 


^(T-Z)  = 


cr+1 


=  b. 


Now,  putting 

■Vy  (r+T -  Z)  =  , 

(nil  =.  b  ,  _ 

2n/y  ^  2*^/^  '■  ‘ 

in  Equation  (5.  5),  taking  the  real  parts,  and  transferring  Z^  to  the  left- 
hand  side,  one  derives  an  explicit  expression  for  Z^  as  follows: 


A 

Z,(Z,  T)  =  —  - 

i'  ’  '  8  n 


(cr-1)^  /( 

^  I /erf(a2) -erf(a^)  +  ^ 


n/tT 


.L\  -y^\ 

i  '  'sin  Za^Y 


-i4y^) . 

-e  siriiiaj^y 


r  b, 


dy\  sin  T  +  / 


n/jt 


Z  Z\ 

- 1^1  -y 


cos  2ajy 


A  1 

1 

(cr+l)2  /[ 

2  'y  , 

'  cos  2a^y 

dy 

>  cos  T 

bn 

J 

J 

n/JT 


e  '  'sinZa^y  -e  '  sinZa^y 


erf(a2,)-erf(a^) 


dy)  sinT 


-5Z-. 


/ 


p  '  fr'rxG  \7-  r1\r  _•»  '  /-. 


n/tT 


'cos  2^27  dy  -e  '  ^  'cos  2ajy  dy 


cos 


1  +  ^  N^r 


\ 


+  -^ 
8 


'  'cos  2a^b^  -  e  ' 

-(4-4)  '  -f* 


COS  Za^b^ 


sin  T 


2  ,2' 

"^1 1 

e  ‘  'sin 2a^b^  -  e 


sin  Za^b 


cos  T(  e 


(^-1) 


^2  ”^2 

cos  Za^b^  -  e  '  ^cos  Za.^h^ 


sin  T 


-W 


e  '  'sinZa^b^  -  e 


H-' 


“1 

\ 

1 

,  (<r+i)A 

'si’n  Za^h^ 

•COS  TJ 

4-^ 

If 


n 


2 


-  a 


m 


./) 


sinZa  y  dy  =  F(a  ,  b  ) 
m'  '  '  m’  n' 


r  -  (4-/) 


^ItT 


cos  2a  y  dy  =  f(a  ,  b  ) 
m'  '  '  m’  n 


then 


Zj(Z,  T)  = 


8  NY 


(cr-1)' 


|j'erf{a2)-erf(a^)+F(a2,  bj)-F(a^,  bj)j 
+  [f(apb^)  -  fta^bpj  cos 


+  e 


(i+ill 

/[erf(a2)-erf(ap+F(a2,  b2)-F(a^  ,  b^)!  sin  T 


T 


J 


sin  T 
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'  ^  /r-faf-bA  -4^-bf) 

1  +  —  fv/JT  \  e  V  'cos  2a  b  -  e  '  ^ 

^  8  X:  11 

r-(a^b^)  -fa^b^) 

+  [e  W  VsinZa^bj^  -  e  ^  ?.a 


cos  Za^bJ  sin  T 


A  [  -  K-»2 

+  -g-  \fir\  e  '  'cos  2a^b^  -  e 


1  )  - 

sinZa^b^  cos  T  /  e  ^ 


cos  Za^b^J  sin  T 


f  ■  ■  ■f^^2)  . 

-  I^e  '  sinZa^b^  -  e  '  'si 


1  i 

sin  Za^b^J  cos  T  f  e  '' 


2\-l 


(5.6) 


Now  putting 


(r  +  T-Z+Zp  =  a, 
i^Ff (T  -  Z  +  Zj)  =  a^ 


-  42  .b^) 

\  .m  nj  . 

i  '  'SI 

-  42  -b^) 

V  m  n/ 

»  '  f  nr\ 


sin  Za  b  =  G(a  ,  b  ) 
rn  n  m'  n 


cos  2a  b  =  g(a  ,  b  ) 
m  n  m  n' 


erf(a  .)-erfa  =  E(a  -a  ) 
'  m'  n  '  m  n' 


we  obtain 


az,  .  , ,2 

- -(Z,  T)  =  —  -  e"^  1 

dT  8  *  Y 


1  ^<^^/7[g(a4,  bi)-g(a3,  bi)j  +f(a3,  bi)-f(a^,  bi)|  sinT 
+  ^'^|G(a^,  bj)-G{a3,  bj)J  -E(a3-a^)  +  F(a^,  b^)  -  F(a3,  b^)^  cos'll 

-e  ^  f<^[g{a3,  b2)-g(a^,  b^)]  +  f(a^,  b^) -f(a3,  b2^  sin  T 


sin  T 
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1 

8 


+  e 


a^,  h2)-G(3.y  b^)]  +  E(a2-a^)+F(a2,  b2)-F(a^,  b^)^  cos  jj 
.  bj)-g(a^,  b^)j  sin  T  +  [G(a2,  bj)-G(a^,  b^)j  cos  t| 
gla^,  b2)-g(a^,  b^)]  smT+  [^G{a^,  b^) -Gia^,  b^)] 


A-1 


cos  r, 


(5.  7) 

In  Equation  (5.  7),  the  term  occurs  on  the  right-hand  side.  To  obtain 
9Z^/8T  at  a  particular  Z  and  T,  therefore,  the  value  of  Z^  has  to  be 
obtained  from  Equation  (5.  6)  for  the  given  Z,  T,  and  then  has  to  be  sub¬ 
stituted  in  Equation  (5.7).  Thus,  the  ratio  c'^'^o  calculated. 

Now,  from  Equation  (5.5) 


A  ft/  '^1 

T7  =  -8  Jv 


-b 


.  t>j^)-g(a^,  bpj  sin  T  +  [^-G(a^,  bj)+G(a2,  b^J  cos  ^  n/7 


e  ^  ^|'g(a^,  b2)-g(a2,  b^)!  sin  T  +  [^G(a2,  b^) -G(a^,  b^)]  cos 'll 

/ 

3.  bj)-g(a^,  bp]  sinT+  jG(a^,  bp -G(a2,  bp]  cos  n/T" 


-1 


(5.8) 


Since 


9Z  9Z 

TZ  ^ 


dT 


9Z^ 

8Z 
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then 


dT 


/  b 

F  =  ^|f(a3.bj)-f(a4.bj)]  sin  T  +  [-E(a3.a4)+F(a4.bj) -F(a3.bp] 


cos 


-bn 


-e 


^f(a^,b2)-f(a3,b2)J  sin  T  +  [E(a3-a^)+F{a3,b2) -F(a^.,b2)J  cos 

(5 . 9) 

The  relative  simplicity  in  the  expression  for 

az. 


JL  ^  1 

u  ax 

o 


is  Observed. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  ballistic  theory  of  an  electron  beam  has  been  developed,  for 
multiple  signals,  subject  to  approximations  made  to  simpfly  solution. 

Some  work,  using  nonlinear  space-charge  wave  analysis,  has  also,  been 
done  and  a  comparison  of  the  ballistic  and  space -charge  wave  analysis 
should  be  attempted  to  throw  more  light  on  the  problem.  The  graphic  plot 
obtained  (Figure  9)  which  gives  a  value  of  G.5  for  the  bunching  parameter 
indicates  that  there  is  pulse  distortion.  The  saturation  of  the  lower  half 
of  the  envelope  indicates  that  the  exit  current  is  rich  in  harmonics,  even 
for  low  depths  of  modulation. 

The  theoretical  study  indicates  that  experiments  can  provide  solu¬ 
tions  where  theoretical  formulation  would  be  cumbersome.  An  experimental 
verification  of  the  theory  should  therefore  be  undertaken,  using  a  set  of 
parameters  designed  to  approximate  closely  the  assumption  of  an  infinite 
beam  and  no  space -charge  effects. 

Attention  is  now  being  directed  to  the  generation  of  high  peak -power 
radar  using  nanosecond  pulses  of  the  Gaussian  type.  Once  such  pulses  are 
generated,  the  response  of  the  klystron  to  them  can  be  observed.  Good 
microwave  amplification  of  these  pulses  would  lead  to  their  use  in  radars. 
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APPENDIX  A:  A. NOTE  ON  THE  GAUSSIAN  SPECTRUM 


As  this  report  deals  with  the  klystron  response  to  a  Gaussian  enve¬ 


lope,  a  short  note  on  the  Gaussian  spectrum  is  useful..  The  analysis  pre- 

at^ 

sented  has  been,  mostly  based  on  the  envelope  g(t)  =  Ve  coswt: 

,2  -ja)„t 


a  • 


g(w  )  =  Re  J 

^  °  Zir  J 

-a 


=  Re 


a 

-i-  f 

2ir  J 


^e-at“^  -o"  ^+ja)t 


(t  +J, 


Ve 


dt 


Thus,  in  the  frequency  plane  the  envelope  is  also  Gaussian.  The 
Gaussian  envelope  is  evidently  economical  in  bandwidth  for  a  given  pulse 
length,  the  majority  of  the  energy  being  confined  to  a  finite  range  of  the 
frequency  spectrum  centered  on  the  carrier  frequency.  These  factors  to¬ 
gether  with  the  fact  that  a  Gaussian  pulse  is  easier  to  generate,  are  the  cri¬ 
terion  determining  its  selection  for  analysis.  It  might  be  interesting  to  note 
to  what  extent  the  envelope  sha.pe  depends  on  the  parameter  a.  Consider 
just  the  envelope  given  by  the  equation. 


f(t) 
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then 


g(<-o) 


!_ 

2 


'v/ira 


Select  a  =  ir/a. ,  and  multiply  both  f(t)  and  by  l/'\[a  ;  then 


TTt 

a 


fl(t)  = 


\fa 


aoi 


gl(«o)  = 


4tt 


Ztt 


The-area  under  the  curve  fj(t)  =  e~  ^  J  va  is  given  by 

,ii 


A  =  2ir 


Zv 


w  =0 
o 


2tt 


Therefore  as  a  becomes  smaller  and  smaller,  the  curve  fj^(t)  becomes 
taller  and  narrower  and  approaches  a  unit  impulse  as  a  approaches  zero. 
Since  a  is  inversely  proportional  to  a,  one  must  have  a  high  value  of  a 
to  obtain  short  pulses.  If  the  frequency  is  increased,  the  value  of  a  has 
to  be  increased  also,  if  the  pulse  is  to  decay  to  a  fixed  fraction  of  its  ampli¬ 
tude  after  a  fixed  number  of  r-f  cycles.  Actually  for  this  purpose  the 
ratio  a/w^  ,  where  u  =  2irx  frequency,  has  to  be  maintained  constant. 

The  transition  from  a  frequency  spectrum  consisting  of  a  series  of 


discrete  frequencies  to  one  consisting  of  a  continuous  band  of  frequencies 
can  be  made  by  treating  the  nonperiodic  function  as  a  periodic  function  in 

at2 

which  the  period  approaches  a  .  The  unit  Gaussian  envelope  e  cosut 
will  be  considered.  The  amplitude  of  the  spectrum  at  ,  is 


g(«o) 


=  A  ..  f 

Ztt  j 

-a 


2.  t 

e-at 


For  a  single,  puls  q,  W:here  .  f(,1;)  =  Q,  for  aU  yalu,es  ,pf  t  except  -L<t<L, 
we  have 


(w  -  u  ) ' 
o' 


L 


g(«o)  =  e  Re  /  e 


t- j 


(«  -  w^) 


i2 


2a 


dt 


^  e  Re 


2ir  ^/4 


/  (w-w  )\  / 

erf^ya  (L  -  j  v  ]  +  eri^l&  (L  +  j 


2a 


where  vise  has  been  made  of  the  identity, 

eff(-x)  =  -erf(x) 

With  the  expansion  for  the  complex  error  function  and  then  taking  the  real 
part,  one  obtains 


g(w  )  =  e  (2  erfN/a  L) 

°  Z  -tt  sla 


ferfi^  g  4a 

\  Tt^la 
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The  -spectrum  of  a  train  of  Gaussian  pulses  of  length  2(L  +  AL)  recurring 
every  T  seconds  -will  be  foun.d  from  the  spectrum  of  a  single  pulse  of  the. 


train.  For  the  single  pulse  at  any  frequency  w^/2tt, 


(to  -  to^) 

.  ^  erf  n/T (L  +  AL)  ^  4S 

°  itn/S- 


For  a  period  of  such  pulses  recurring  with  a  spacing  T  =  l/C  ,  the  sum 
of  spectra  of  the  individual  pulses  form  a  Fourier  series  of  harmonics  of 
C  ;  therefore 


f(t) 


A  + 
o 


cos  2TTnCT 


where  is  the  sum  of  an  infinite  number  (one  from  each  pulse)  of  infini¬ 

tesimal  terms  g(2TrnC)  and  g(  -2TrnC)  ,  giving 


A 

n 


^ — •  (to' -  2iTnC)^ 

\  erf 'vya  (L  +  AL)  4a 

Tt'sla 


To  put  an  absolute  value  on  the  amplitudes  >  it  is  necessary  to  aver¬ 

age  them  over  the  recurrence  period  of  the  single  pulse,  making  them  in¬ 
finitesimals,  However,  in  the  train  of  pulses  recurring  every  T  =  l/C 
seconds,  the  amplitude  of  A^  can  be  determined  by  averaging  the  terms 
in  g(w^)  over  an  interval  T;  then 


A 

n 


erf  N/a  (L  +  AL) 


(to  -  2TmC) 
4a 


■ir  ^  T 


e 


and  when  T  =  4L  + 


1 

C  ' 


thus, 


where 


C  erfN^  L  fl  + 

A  =  _ ^ _ i _ L 

Tr'vSr(4LC +■•  1, 


(w  -  ZirnC)^ 
4a 


f(t)  =  ’*'■^1  cos  ZirCt  +  A^  cosZirZCt  +  .  ,  . 


A  is  known, 
n 
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APPENDIX  B.  FOURIER.  COEFFICIENTS 


The  Fourier  coefficients  and  b^, ,  as  obtained  in  Equation  (2,  9) 


are  simplified  here: 


a  =  —  /  I  cos r 

r  IT  j  o 


y  -  k  sin  my  + - -  k  sin  ny  +  — 


dy 


It  is  sufficient  to  obtain  the  solutions  to  the  following  coefficients,  as  they 
are  related  to  the  ones  in  Equation  (2.  9): 


IT 


=r  y  cos  r  sin  (my  +  +  c^  sin(ny  +  a^)!  dy  , 


IT 


a-  +  j 

r  r 


a'  -  j  b' 
r  r 


^  / 


y  sin  r  j'y  +  sin  (my  +  a^)  +  sin(ny  +  a^)!  dy  , 
■ir 

^  jr  ^y  sin  (my  +  +  c^  sin(ny  +  dy 


=  I 


r  -j  r  ^y  +  Cj  sin  (my  +  +  c^  sin(ny  + 


Now,  according  to  the  property  of  Bessel  Functions, 


ejzsine  ^  +  2 

^2^^)  20 

+  J^{z) 

cos  40  +  ..  .. 

+  2  j 

^J^(z)  sin  0 

+  J3(z) 

sin  30  +  .  .  , 

g 

1 

=  jp(z) 

ejP® 

(B.  1) 

p=-a 
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since  J  -  p(z) 


Jp(z)  ;  therefore, 


using  the  property  in  Equation  (B,  1) 


gives 


I 

a 

r 


+  Jb, 


a 


jp(my +  a^) 
e‘ 


(rc^) 


j  q(ny  H-a^) 
e 


dy 


p=-a 


q=-a 


p=-a  q=- 


Jp{rCi)  Jqirc^) 
a 


j(pm  +  qn)  y  +  jCpUj  +  qa2) 
e  dy 


V 


Because  of  the  nature  of  the  integrand,  the  order  of  integration  and 
summation  can  be  interchanged;  therefore 
a  a 

j{p°^2+ sin(pm  +  qn  +  r) 


j*^r  =  ^  J 


Jp(rci)  J  (fc^)  e 


pm  +  qn  +  r 


p=-a  q=-a 


Now,  noting  that 


J  (-z)  =  J  (ze'^^)  = 

V  '  V  ' 


t  ik 

_  m=o _ ^ 


.  \v+Zm 
zeJ"^) 


-  >  \2 


v+2m 


mi  r  (v  +  m  +  1) 


mj  r  (v  +  m  +  1) 


m=o 


j  Zirm  .  - 

since  =  1  ,  therefore 


J  (-Z)  =  J  (z) 

V  v'  ' 


Using  the  same  procedure  as  before,  we  obtain 
a  CL 

,  ^ 
a 

.  M  J. 

p=-a  q=:-a 


T 


S5 


r 
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APPENDIX  C.  EVALUATION  OF  erf{a  -  jb) 


In  the  analysis  of  klystron  response  to  Gaussian  wave  excitation, 
cornplex  error  functions  of  the  type  erf(a  -  jb)  ,  where'  a  and  b  are  real 
numbers,  have  often  been  encountered.  Here,  a  method  will  be  indicated 
to  express  the  complex  error  function  in  terms  of  real  integrals.  Consider 
the  expression, 

a-jb  . 


erf(a-jb)  = 


N/ir 


/ 


e"^  dz 


f(z)  dz  =  0  around  a  closed 


in  the  complex  plane  with  z  =  x  +  jy  . 

2  f 

As  e"2  is  an  analytic  function,  ® 

path  of  integration.  Hence,  integrating  along  the  path  shown  by  the  arrows 

in- Figure  13  gives 


a  p  D 

j  e'^  dx  -  j  J 


e-(a-jy)  dy 


a-jb 


-I 


-z 


dz 
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Figure  13,  Path  of  Integration  for  the  Function  e” 
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from  which 


X 


\ 


> 


i 


1 


erf{a  -  jb)  = 


r  ~x^  .  2  1  -{a^  -  y‘')  +  2jay 


2  2, 


—  I  e  dx - ;  e 

n/tt  J  n/iT 

o  o 


dy 


b  b 

=  erf(a)  +  — ^  /  ^  sin2ay  dy  -  j  — ^  ^  cos  2ay  dy 

/  'v/if 


This  representation  of  a  complex  error  function  in  terms  of  real  integrals 
has  been  used  frequently. 

Some  other  important  results  follow: 


erf(-a-jb)  =  erf(-a) 


U 

— —  I  e’^^  ^  sinZay  dy  -  ■—  ^  cos  2ay  dy  , 

\flr  J 


as 


erf(a)  =  -erf(-a) 


a-jb 


erf(a-jb)  -  erf(-a-jb)  = 


2 

's/ir 


e  dz 


-a-jb 


=  2  erf(a)  +  —  /  e”^^  "  V  )  sin  2ay  dy  , 


which  is  a  real  quantity.  This  same  result  can  be  obtained  by  contour 
integration  around  a  suitable  rectangle  in  the  z -plane.  Again 
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erf(a  -  jb)  -  erf(c-jb) 


=  erf{a)  -  erf(c)  + 


\fn 


b 

/ 


[e-(^ 


y^) 


sin  2ay 


These'  results  show  that  a  complex  error  function  can  be  easily 
computed,  and  its  real  and  imaginary  parts  separated.  Considerable  sim¬ 
plification  in  computation  can  result  in  specific  problems.  The  asymptotic 
expansion,  for  example. 


erf{a) 


(Za^) 


1.-3 
2  '2 


I- 3,- 5 

2  ^ 
(Za'^) 


is  convenient  for  computation  when  a  is  large. 
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